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Field of the InverLfcior. 

This l^ntion relate, to a ™thod of processing .narine 
-aonetio »radi«xt data and also to «^loratlon methods 
usxng that data. 

Backgroxind Art 

Exploration of the earth to locate oil, minerals or other 
valuable deposits, uses seismic techniques in order to 
provide information about the subsurface structure of the 
earth so that an assessment can be made as to whether a 
particular deposit may or may not be present. m 
particular, seismic techniques provide an indication of 
various subsurface structures including porous sandstone 
or fractured carbonates which may contain oil or other 
valuable deposits. 

20 Whilst this technique does provide an indication of 
subsurface structures, it is possible a particular 
structure which may appear to be of interest, in fact 
represents a significant hazard to drilling. For exa^le, 
xt xs possible that subsurface volcanos e^ibit very 

25 sxmilar characteristics to anticlines which may contain 

oxl. If a drilling operation is initiated and it is found 
that the anticline is in fact a volcano, then a 
considerable financial loss results because of the cost of 
establishing the drilling platform. This is particularly 
3 0 the case in relation to marine e^loration because of the 
sxgnxficantly higher drilling cost. 

in order to provide more information concerning the nature 
Of subsurface structures, magnetic data of the e^loration 
35 area is obtained. 
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A magnetic gradient survey enables magnetic data relating 
to the survey region to be obtained and this magnetic data 
can be used to provide information relating to the nature 
of subsurface structures. if the magnetic survey is 
overlayed with the seismic survey, then structures which 
appear to be of interest from the point of view of the 
seismic survey, can be further considered in the light of 
the magnetic data so that a clearer indication can be 
formed as to whether the structure is a structure which 
may contain a deposit of interest such as an oil deposit, 
or whether the structure exhibits magnetic phenomena, such 
would be the case with a volcano. Thus, the location of 
drilling platforms can be decided with more precision to 
avoid subsurface structures which could be hazardous from 
the point of view of a drilling operation. However, 
conventional techniques for processing magnetic gradient 
data contain considerable distortion and anomalies due to 
unwanted magnetic effects including the ship bias referred 
to above. 

oaxe conventional method of obtaining that data in relation 
to onshore e^loration is to tow magnetometers behind an 
aircraft to obtain magnetic data in relation to the survey 
area. in the case of marine exploration, the 
magnetometers are towed behind a ship. 

The usefulness of marine magnetic data which is obtained 
xn this manner has been limited by the quality of the data 
which is obtained. One of the major problems with marine 
magnetic data acquisition is interference or so-called 
ship bias which is created by the magnetic field induced 
by the ship which tows the magnetometers. 

in conventional marine magnetic data acquisition, two 
magnetic field sensors, which are generally referred to as 
fish, are towed behind a ship. The magnetometers are 
connected to a tow line and the magnetometer which is 

^^«.l..tu^^^su^.^^,^,,^ - Ke«»a oe ^^^^ 
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closest to the ship is towed at a distance of some 300-600 
metres behind the ship to avoid the effect of ship induced 
magnetic field. Furthermore, the magnetometers are 
separated by a distance of in excess of 100 metres. The 
reason for the length of the tow line and the separation 
of the magnetometers is to reduce ship bias and therefore 
provide data which is relatively free of that bias. 

However, because the towing line is so long and the 
distance of separation so great, the degree of drift of 
the magnetometers in the ocean as the magnetometers are 
towed behind the ship is considerable. Furthermore, the 
processing techniques used to acquire magnetic gradient 
data assume that after a given time period, the trailing 
magnetometer will be towed to a position which coincides 
wxth a previous position of the magnetometer which is 
closest to the ship. This assumption is made during the 
processing of the magnetic data. However, because the tow 
line is extremely long and the distance of separation 
between the magnetometers is considerable, drift of the 
magnetometers due to sea currents and the like means it is 
unlikely that the trailing magnetometer will actually 
occupy the same position as the first magnetometer after 
that given time period. 

Furthermore, in the acquisition of the magnetic data, the 
ship is required to travel along predetermined survey 
Ixnes and because of the drifting of the sensors, the 
sensors do not actually travel along the survey line. m 
the southern hemisphere, if the magnetometers drift to the 
north of the survey line, the recorded gradient data 
between the two magnetometers will have larger distortion 
than if drifting towards the south of the survey line. 

35 Summary of the Invention 

The object of the invention is to provide an improved 
processing method which accounts for ship-induced bias and 

vmr mv . itethod of Preeasslng ttoriaa KBoiMtic Oradtane Itata.dee 21/12/01 
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which therefore enables the magnetometers to be towed by a 
tow line which is nmch shorter than the conventional tow 
line and also be located closer together to reduce the 
amount of drift and also increase the likelihood that the 
trailing magnetometer will occupy the same position which 
was previously occupied by the leading magnetometer after 
a given time period. 

The invention may be said to reside in a method of 
processing marine magnetic data obtained by towing first 
and second spaced apart sensors behind a ship, the method 
including; 

obtaining raw magnetic gradient data from the 

sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors/ 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
20 a data output. 

Because the method determines the trend of the gradient of 
the ship bias, a more accurate estimation of the actual 
gradient bias introduced by the ship, having regard to any 
25 drifting of the sensors, is taken into account. By 

subtracting the trend of the gradient bias from the raw 
gradient data, the bias introduced by the ship can be 
removed from the data so that more accurate data is 
obtained. Because the present method eliminates the 
effect of the ship bias, the sensors can be towed much 
closer to the ship and at much closer separation 
distances, thereby reducing the amount of drift of the 
sensors and therefore resulting in the sensors being more 
Ixkely to follow the actual survey line upon which the 
ship travels and also for the rear sensor to occupy the 
same position at a later time which was previously 
occupied by the front sensor. Thus, the ^ality of the 
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output data is greatly improved and therefore provides 
data which gives a much more accurate indication of the 
subsurface magnetic characteristics of the survey area, 
which can be then used in combination with seismic data to 
make an assessment as to the viability of further 
exploration or a drilling program. 

Preferably the method includes the step of determining an 
estimate of the gradient of the ship bias from the raw 
magnetic gradient data obtained by the sensors, and 
determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 

Preferably the raw magnetic gradient data is calculated as 
the measured magnetic signal at the leading sensor less 
the measured magnetic signal at the trailing sensor, 
divided by the distance between the sensors. 

The estimate of the gradient of the ship bias may in fact 
siii«)ly be derived from the raw magnetic gradient data 
which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects. 

Preferably the trend of the gradient of the ship bias is 
determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8. 

Preferably the processing coa^rises the application of the 
corrected gradient data over sample intervals, integration 
of the corrected gradient data into total magnetic 
intensity data, and application of a low pass filter to 
the integrated total magnetic intensity data to provide 
the data output. 

Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 



the following elation: 



TMI=|; [GcAx,(f)]+M,«„ 

i=l 

Wherein Go is the corrected gradient data 
obtained after subtraction of the trend of the gradient of 
the ship bias from the raw gradient data, and 

Mstat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
other . 

The values of TMI so calculated are then smoothed to 
remove residual noise using a low pass filter. 

Preferably the gradient of the ship bias is calculated 
according to the following equation: 

Gm.s = ^/(^)-^^(^) _ ) - Dif^ )] + [M, (f, ) - (rj] 

where : 

M/x) = Meix) + Z)(/i) + M*(ri) 
MXjc) = AfeCx) + D{fii + Mbifii 

where 1«£ is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
Me(x), diurnal variation D(ti), ship bias I*,(tx) caused by 
ship-induced field, sensor drifting and sensor 
perturbation, all at time ti and along track distance x, 
and at some later time ta, Al is the distance between the 
sensors, D(t2) is the diurnal variations sensed by tbe rear 
sensor, and 

Mb(t2) is the ship bias at time ta detected by the 
rear sensor 

\\Mll»_eil«s\hcna$\biiisa\K«>p\Sp«:lXBMP Prov - Method oC Freeeseliig Kariaa Hasaatie CrwUent Itata.doc 21/13/01 
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Ite invention may also b« said to reside in a .athod of 

-» exploration proara», tj» 

method including: 

towing first and second sensors behind a ship 
alonsr predetermined survey lines; 

obtaining raw magnetic gradient data from the 

sensors ; 

determining the trend of the gradient of the ship 
bxas detected by the sensors; 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
a aata output. 

Preferably the method includes the step of determining an 
estxmate of the gradient of the ship bias from the raw 
^gnetic gradient data obtained by the sensors, and 
determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 

preferably the raw magnetic gradient data is calculated as 
the measured magnetic signal at the leading sensor less 
the measured magnetic signal at the trailing sensor, 
divided by the distance between the sensors. 

^ estimate of the gradient of the ship bias may in fact 
sxmply be derxved from the raw magnetic gradient data 
whxch xs Obtained from the sensors and which may include 
data relatxng to the gradient of the ship bias as well as 
data relatxng to other magnetic effects. 

Preferably the trend of the gradient of the ship bias is 
determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8. 

Preferably the further processing comprises the- 
application of the corrected gradient data over sa^le 
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intearvals, integration of the corrected gradient data into 
total magnetic intensity data, and application of a low 
passed filter to the integrated total magnetic intensity 
data to provide the data output. 

Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 
the following eqpiatxon: 



10 TMI=£ [GcAx,it)]+Ms,at 

1=1 

wherein Gc is the corrected gradient data 
obtained after subtraction of the ship bias trend from the 
raw gradient data, and 

Matat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
other. 



15 



The values of OMZ so calculated are then smoothed to 
20 remove residual noise using a low pass filter. 

Preferably the gradient of the ship bias is calculated 
according to the following ecjuations 

25 GMas = ^/<^>-^^(^)_ [g(f,)-P(f,)] + [M,(r,)-Af,(r,)] 

A/ A/ 



30 



where : 



M/x) = Me(x) + D(ti) + Mbih) 

MAx) = Me(x) + £)(/2) + Affr(?2) 



where Mf is the magnetic field as measured by the 
front sensor and which consists of the envirozimental field 
35 Me(x), diurnal variation D{ti), ship bias Mb(ti> caused by 
ship-induced field, sensor drifting and sensor 

XX«.U._Ml«.Vho»eS\U.l«.XK«pVSp«lMWP p™, . of Pr^Maln, «a,i„« M«,n-tlc «r«U«t Oacdoc ai/ia/oi 
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perturbation, all at time ti and along track distance 
and at some later time tj, Al is the distance between the 
sensors, D(t2) is the diurnal variations sensed by the rear 
sensor, and 

Mb(t2) is the ship bias at time ta detected by the 
rear sensor Mr* 

The invention may also be said to reside in a method of 
subsurface exploration to determine viability of drilling 
in a marine environment by considering magnetic data 
relating to the environment, and wherein the magnetic data 
has been obtained by: 

obtaining raw magnetic gradient data from the 

sensors; 

15 determining the trend of the gradient of the ship 

bias detected by the sensors; 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; 

processing the corrected gradient data to provide 
20 a data output. 

Preferably obtaining the magnetic data includes the step 
of determining the gradient of the ship bias from data 
obtained by the sensors, and determining the trend of the 
25 gradient of the ship bias from the gradient of the ship 
bias. 

Preferably the raw magnetic gradient data is calculated as 
the measured magnetic signal at the leading sensor less 
the measured magnetic signal at the trailing sensor, 
divided by the distance between the sensors. 



30 



The estimate of the gradient of the ship bias may in fact 
simply be derived from the raw magnetic gradient data 
35 which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects. 

\Mi«ib.£ll«\ho«e$\Uil«i\Keep\Speel\BHP Prov . Method ot Proeesaing Itorino HiMHiecie Gradient Oate.doc 21/12/01 
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Preferably tlie trend of the gradient of the ship bias is 
determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8. 

5 

Preferably the processing comprises the application of the 
corrected gradient data over sample interval S/ integration 
of the corrected gradient data into total magnetic 
intensity data, and application of a low passed filter to 
10 the integrated total magnetic intensity data to provide 
the data output. 

Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 
15 the following equation: 

wherein Gc is the corrected gradient data 
obtained after subtraction of the ship bias trend from the 
20 raw gradient data, cuid 

Mstat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
other. 

25 

The values of TMX so calculated are then smoothed to 
remove residual noise using a low pass filter. 

Preferably the gradient of the ship bias is calculated 
30 according to the following equation: 

Al A/ 

where : 

35 
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M/x) = Meix) + DitO + Mb(ti) 

Mr(x) = Meix) + D(/2) + Mi^Cfj) 

5 where Mf is tlie magnet: Ic £ield as measured by the 

front sensor and which consists of the environmental field 
Me(x), diurnal variation D(ti), ship bias Mb(ti) caused by 
ship-induced field, sensor drifting and sensor 
perturbation/ all at time ti and along track distance 
10 and at some later time t2, Al is the distance between the 

sensors, D(t2} is the diurnal variations sensed by the rear 
sensor, and 

Mb(t2) is the ship bias at time ta detected by the 
rear sensor M^* 

15 

The invention may still further be* said to reside in a 
method of drilling for a deposit in a marine envirozment, 
including: 

determining the location of drilling from data 
20 which has been obtained and which indicates the possible 
existence of the deposit; and 

which location is also determined by magnetic 
data which has been obtained by towing magnetic sensors 
behind a ship, the magnetic data being processed by: 
25 obtaining raw magnetic gradient data from the 

sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 
subtracting the trend from the raw magnetic 
30 gradient data to obtain corrected gradient data; 

and 

processing the corrected gradient data to provide 
a data output. 

35 Preferably the 'method includes the step of determining an 
estimate of the gradient of the ship bias from the raw 
magnetic gradient data obtained by the sensors, and 
determiziing the trend of the gradient of the ship bias 



N\nelb_£iles\ho»e$\LuU«\Keep\Speci\8l» Prov - Method of Processing Katrine Hagnetie Gradient Data. doe 21/12/01 
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from t:liat; estimate of tlie gradient of the ship bias. 

Preferably the raw magnetic gradient data is calculated as 
the measured magnetic signal at the leading sensor less 
5 the measured magnetic signal at the trailing sensor, 
divided by the distcuace between the sensors. 

The estimate of the gradient of the ship bias may in fact 
simply be derived from the raw magnetic gradient data 
10 which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects. 

Preferably the trend of the gradient of the ship bias is 
15 determined by applying the gradient of the ship bias data 
to a B-Spllne filter with smoothness 0.8. 

Preferably the processing of the corrected gradient data 
comprises the application of the corrected gradient data 
20 over sample intervals. Integration of the corrected 

gradient data into total magnetic intensity data, and 
application of a low passed filter to the integrated total 
magnetic intensity data to provide the data output. 

25 Preferably the total magnetic intensity is obtained by 

Integrating the corrected gradient data in accordance with 
the following equation: 

30 wherein Gc Is the corrected gradient data 

obtained after soabtractlon of the ship bias trend from the 

raw gradient data, and 

Mstat is the total magnetic intensity value at a 

point representing the start of the survey or at a 
35 location where a survey line and a tie line cross each 

other . 

\\melb.£ilos\boae$\Uaisa\Ke«p\5peel\BRP Prov - Method of PceeesAing Karine Hagnaeie Gcadleae Daea.dee 21/12/01 
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The values of TMI so calculated are then smoothed to 
remove residual noise usingr a low pass filter. 

Preferably the gradient of the ship bias is calculated 
according to the following equation: 

Gi,i„ = ^/(^)-^rW CZ)(/,)-D(r,)] + [M,(^,)-M,(f,)] 

Ai 

where : 

Mj(x) = Meix) + D(ti) + Mbiti) 

Mr(x) = Me(x) + D(t2) + M^itz) 

where Mf is the magnetic field as measured by the 
front sensor and which consists of the envirozimental field 
Ma(x)^ diurnal variation D(ti)^ ship bias Ife(ti) caused by 
ship-induced field, sensor drifting and sensor 
perturbation, all at time tx and along track distance x, 
and at some later time ta, Al is the distance between the 
sensors, DCta) is the diurnal variations sensed by the rear 
sensor, and 

Mb(t2) is the ship bias at time t2 detected by the 
rear sensor Mr. 

Brief Description of the Drawings 

A preferred embodiment of the invention will be described, 
by way of example, with reference to the accoxnpanying 
drawings in which: 

Figure 1 is a view showing a marine magnetic data 

survey/ 

Figure 2 is a flow chart showing the preferred 
method of processing the marine magnetic gradient data 
according to the preferred embodiment of the invention; 

Figure 2A is a graph showing the function of a 
low pass filter used in the preferred embodiment; 

\\«aelb.eHo«\ho«e$\Lttio«\ltoep\SpeciNaHP Prov - Method of Procesoiag Marine tfagnecie Cr«U«.e Data.doc 21/12/01 
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Figtixe 3A and Figure 3B are graphs showing data 
obtained according to the survey and as processed 
according to the preferred embodiment of the invention; 

Figure 4 is a representation of a survey using 
5 data created according to a prior art technique; and 

Figure 5 is a diagram similar to Figure 4, but 
using data obtained according to the preferred embodiment 
of the present invention. 

10 Description of the Preferred Embodiment 

With reference to Figure 1^ a marine survey to acquire 
magnetic data is illustrated. A ship 10 tows first 
magnetometer Ml smd second magnetometer M2 by a tow line 
12. a*he ship travels along survey line 14 and the line 

15 16/ which are dotted, represent the ship- induced magnetic 
field. 

In Figure 1, it is assumed the ship is in the southern 
hemisphere and is heading north-northwest, and the sensors 
20 are drifting towards the northeast. In this case, the 

recorded gradient data between sensors Ml and M2 will have 
a larger distortion than if the sensors were drifting 
towards the southwest. 

25 Ideal conditions would assume that the sensors Ml and M2 
are towed directly behind the ship on the survey line auxd 
that therefore, after a given time period, the sensor M2 
would occupy the same position as was previously occupied 
by the sensor Ml. However, due to drifting of the sensors 

30 because of ocean currents and the like^ this ideal 

situation in general cannot be achieved. Nevertheless, 
the processing methods used to process magnetic data do 
assume that the sensor M2 will occupy the same position as 
the sensor Ml, previously occupied, but at a later time. 

35 The processdLng method according to the preferred 

embodiment en ables ship bias to be accounted for much more 
accurately than any prior art techniques axid therefore 

\\molb.fil«s\homo$\Lttiu\Ke«p\Specl\BHP Pm - Hdchod of Frecesaing Marine Hagneeic Gradieae Oata.dee 21/12/01 
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enables the sensors to be towed closer- to the ship and to 
be separated by a smaller distance. This therefore 
provides much greater accuracy in the measurements because 
the amotint of drift is not as great due to the shorter tow 
line and also because the separation between the sensors 
is smaller, the likelihood of the rear sensor occupying 
the same position as the front sensor at a later period in 
time is greatly increased. 

At any time ti and along-track distance x, the front sensor 
measures magnetic field - Mf consisting of the environment 
field Uelx), the diurnal vsuriation D(ti), ship bias Jkfe(ti) 
caused by ship induced field, sensor drifting and sensor 
perturbat ion . 

M,(x) = Meix) + Z>(f i) + Mbiti) 

At some later time tj, the rear sensor makes a measxur^aent 
at the same spatial point: 

M^X) = MeOO + D{t2) + Mb{t2) 

The difference between the two measurements at the same 
location is: 

M/Oc) - MA.X) = [Ditx) - D(t2)] + [Mb(ti) - Mb(t2)] 
Ship bias gradient (Step 2, Figure 2) can be expressed as: 

A/ A/ 

It should be mentioned that the calculated ship bias 
gradient is still affected by diurnal variation. In 
practice we notice that the sea current drifting is a kind 
of long wave-length variation so ship bias should be 
coincident with such current drifting. step 2 in Figure 

\\m.ll».«iX«.\homeS\Lia.«Me«.p>Speei\8l» Pcov - Itochod of PrecMing K.rin. M^rMcie Or«Uau: D.t..doe 21/«/01 
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2, which is the calculation of the srradient of the ship 
bias, is actually an estimate of the gradient of the ship 
bias which is obtained from the raw data collected by the 
sensors Ml and M2. The data collected by the sensors will 

5 include raw magnetic gradient data which contains many 

signals including environmental magnetic signals, diurnal 
signals and the ship bias, as well as instrument bias and 
drift. The trend of the gradient of the ship bias is 
obtained from this estimate by applying the estimate of 
10 the gradient of the ship bias to a B-Spline filter, as 

will be described in more detail hereinafter. Therefore 
the trend of the gradient of the ship bias (Gtrend) (Step 

3, Figure 2) is used for the correction along a track 
line. 



15 



20 



25 



35 



The gradient correction of the ship bias {Gt)oaxk be 
expressed as follows: 

Gc = G- Gtrend (Step 4, Figure 2) 
In this equation, G is the raw magnetic gradient data. 



As said before such drifting can cause large bias effects. 
Conseg:uently, when the sensor drifting and perturbation 
take place simultaneously the bias zig-zag will be around 
its mean values. The variable mean value of the bias 
along a survey line is regarded as bias trend since the 
high fre<xuencies of sensor perturbation only create reoxdom 
noise around the bias trend and after integration its 
30 effect aan be eliminated from raw gradient data. 

After removal of ship bias effects TMI (Total Magnetic 
Intensity) can be calculated (Step 6) through integration 
of magnetic gradient data: 



TMI= [GcAxi(t)]+Mstat 



\\Bal^eUas\haina$\bul8a\K«ep\fipeel\BHP Prov - Kaehed at Preesaaiag Hartaa MaoaaCie aradlane 
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Hore^ Axi(t) is tlie sampling distance along the susrvey 
line. Wstat is tlie TMI value at a point of survey start or 
at the location where the survey line and a tie line cross 
5 each other. 

The TMI value calculated in Step 6 is then smoothed using 
a low pass filter (Step 7) to remove any features which 
have a greater rate of change of TMI with distance than is 
10 expected in the particular survey area. An example of the 
action of this filter is provided in Figure 2A in which 
trace 50 is the smooth TMI curve and trace 51 is the TMI 
data prior to smoothing. 

15 The data output obtained at Step 7 may include line 

levelling and data griding to provide a final data output. 

Figures 3A smd 3B represent an actual example of the 
preferred embodiment of the invention which relates to raw 
20 data obtained from a known region* 

In Figure 3A^ the trace 20 represents the total ship bias. 
The line 21 represents the bias trend and it can be seen 
from the left hsoid side pf the trace 20 to the right hand 

25 side of the trace 20^ the bias trend 21 changes non- 

uniformly about the value of -0.08 on the graph of Figure 
3A. A particular value for the bias trend can therefore 
be obtained for various intervals and this value can be 
subtracted from the raw data in order to provide the 

30 corrected data. The trend of the gradient of the ship 

bias is preferably determined by calculating the gradient 
of the bias data 20 and then applying that gradient data 
to a B-Spline filter with smoothness 0.8 in order to 
provide the representation of the ship bias 21. 

35 

In the preferred embodiment of the invention^ the gradient 
of the bias data 20 is calculated simply by using the raw 

\\meU>.eUes\hoBe$\Lul8a\Keep\Speel\aaP Prov - Method of Proceaslng Karino Kagneeie Gradient Data. doe 21/12/01 



- 19 - 



magnetic gradient data obtained from tlie sensors because 
when that data is applied to the B- Spline filter, all that 
remains is the gradient of the ship bias component. The 
ship bias 21 can then be subtracted from the calculated 
5 raw gradient. 

The trend of the gradient of the ship bias is a non-linear 
function which is represented by the line 21 in Figure 3A. 
The line is a measure of how the bias chemges over time as 

10 the sensors are towed behind the ship. As is clearly seen 
from line 21, the trend is not constant or merely an 
average, but rather represents the fluctuation of the bias 
gradient and at some times, is above the value of 0.08 and 
sometimes below that value. As noted above, the trend is 

15 determined by applying the gradient of the ship bias to a 
B-Spline filter with smoothness 0.8. However, in other 
«nbodiments, the filter could have different smoothness 
depending on the region from which data is collected and 
the nature of the data which is collected. In general, 

20 the function of the filter is to smooth the trace 20 so 

that some meaningful value of the bias trend at particular 
time periods can be obtained. Thus, the filter is 
effectively determining the peaks and troughs of the trace 
20 and fitting a curved line between those peaks and 

25 troughs, which gives a measure of how the gradient of the 
ship bias is changing with time as the sensors are towed 
behind th^ ship. 

In the example of data obtained from the known region, 
30 dashed line 30 represents total magnetic intensity data 
which is provided according to the prior art processing 
techni<3ue. The line 32 r^resents the integrated total 
magnetic intensity data from the bias corrected gradient 
data according to the preferred embodiment of the 
35 invention. The line 34 represents the diurnal variation 
at a station some 500 km away from the survey area and 
line 36 is the observed field data including the diurnal 
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effect. 

It can be seen from the conventional processing techni<iue, 
which gives line 30, false anomalies are provided which do 
not show up on the line 32, which is produced according to 
the present invention. 

In general, the results of the magnetic survey are 
produced in a colour chart. Black and white 
representations of a colour cheurt are shown as Figure 4 in 
relation to conventional processing of data obtained from 
a known region and Figure 5 by a processing method 
according to the preferred embodiment of this invention. 

The characteristics of the Igneous basement of the known 
region are well known and it can be seen that data 
processed according to the present invention produces a 
much clearer indication of actual magnetic structures than 
P^^or art technic[ue which include considerable 
interference and data which may mislead emalysers into 
considering that magnetic structures exist which are not 
actually there, or that magnetic structures do not exist 
when a magnetic structure actually is present. 

Thus, the preferred embodiment of the present invention 
provides data which more accurately reflects the 
likelihood of magnetic structures which can then be used 
as a basis for determining the viability of a drilling 
operation in an exploration or recovery program. 

Since modifications within the spirit and scope of the 
invention may readily be effected by persons skilled 
within the art, it is to be understood that this invention 
is not limited to the peurticular embodiment described by 
way of example hereinabove. 
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